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Abstract 

A path from s to t on a polyhedral terrain is descending if the height of a point p 
never increases while we move p along the path from s to i. No efficient algorithm 
is known to find a shortest descending path (SDP) from s to t in a polyhedral 
terrain. We give two approximation algorithms (more precisely, FPTASs) that solve 
the SDP problem on general terrains. Both algorithms are simple, robust and easy 
to implement. 
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1 Introduction 

Finding a shortest path between two points in a geometric domain is one of 
the most fundamental problems in computational geometry. One extensively- 
studied version of the problem is to compute a shortest path on a polyhedral 
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terrain; this has many apphcations in robotics, industrial automation, Geo- 
graphic Information Systems and wire routing. Our paper is about a variant 
of this problem for which no efficient algorithm is known, the Shortest De- 
scending Path (SDP) Problem: given a polyhedral terrain, and points s and t 
on the surface, find a shortest path on the surface from s to t such that, as 
a point travels along the path, its elevation, or z-coordinate, never increases. 
We need to compute a shortest descending path, for example, for laying a 
canal of minimum length from the source of water at the top of a mountain to 
fields for irrigation purpose, and for skiing down a mountain along a shortest 
route [1,17]. 

The SDP problem was introduced by de Berg and van Kreveld [9], who gave a 
polynomial time algorithm to decide existence of a descending path between 
two points. Since then the problem has been studied in different restricted 
settings [1,3,17] (See Section 2.2 for a brief survey), but the SDP problem on 
general terrains remained open in the sense that neither a polynomial time 
algorithm nor a polynomial time approximation scheme (PTAS) was known. 
In this paper we present two approximation algorithms (more precisely, fully 
polynomial time approximation schemes, FPTASs) to find SDPs in general 
terrains. These algorithms have appeared in preliminary forms in Ahmed and 
Lubiw [2] and in Roy et al. [18] respectively. Both the algorithms discretize 
the terrain by adding Steiner points along the edges, thus transforming the 
geometric shortest path problem into a combinatorial shortest path problem 
in a graph. This approach has been used before for related shortest path 
problems such as the Weighted Region Problem and the Shortest Anisotropic 
Path Problem (discussed in Section 2.2). In those results, Steiner points are 
placed independently along each edge. Such independent placement fails for 
SDPs, and our main new ingredient is to place Steiner points by slicing the 
terrain with horizontal planes. Both the algorithms presented here are simple, 
robust and easy to implement. 

In our first algorithm, given a vertex s in a triangulated terrain, and a constant 
e G (0, 1], we discretize the terrain with O (~^) Steiner points so that after 

an O (^^^ log (^)) -time preprocessing phase, we can determine a (1 + e)- 
approximate SDP from s to any point v in 0{n) time if v is either a vertex 
of the terrain or a Steiner point, and in O time otherwise, where n is 

the number of vertices of the terrain, and X is a parameter of the geometry 
of the terrain. More precisely, X = ^ ■ = ^ sec 9, where L is the length 
of the longest edge, h is the smallest distance of a vertex from a non-adjacent 
edge in the same face (i.e. the smallest 2D height of a triangular face), and 
6 is the largest acute angle between a non-level edge and a vertical line. Our 
second algorithm places Steiner points in a different manner, which modifies 
the above preprocessing time and the two query times to O ( ^^^^^ log^ ( ) ) , 



0{n), and O (^^log(^^]) respectively, where X' = ^. In comparison, the 
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first algorithm is faster in terms of n, e and ^, but it depends heavily on the 
inclination of the non-level edges. On the other hand, the second algorithm 
does not depend at all on edge inclinations, and hence is better for terrains with 
almost level edges. It is straightforward to follow a "hybrid" approach that first 
checks the edge inclinations of the input terrain, and then runs whichever of 
these two algorithms ensures a better running time for that particular terrain. 

The paper is organized as follows. In Section 2 we define a few terms, discuss 
the properties of SDPs, and mention related results. Sections 3 and 4 give 
details of our approximation algorithms. We conclude in Section 5 with a few 
open problems. 



2 Preliminaries 

2. 1 Terminology 

A terrain is a 2D surface in 3D space with the property that every vertical line 
intersects it in at most one point [10]. We consider triangulated terrains. For 
any point p in the terrain, h{p) denotes the height of p, i.e., the 2-coordinate 
of p. We assume without loss of generality that all points of the terrain lie 
above the plane z = 0. An edge or face in 3D is level if all points on that edge 
or face have the same height. We add s as a vertex of the terrain. Let n be 
the number of vertices in the terrain. By Euler's formula [10], the terrain has 
at most 3n edges, and at most 2n faces. 

We reserve the terms "edge" and "vertex" for features of the terrain. We 
use the term "segment" to denote a line segment of a path, and "node" to 
denote an endpoint of a segment. We use "node" and "link" to mean the 
corresponding entities in a graph or a tree. Figure 1 shows the convention we 
will use in our figures to mark various components related to a descending 
path. In particular, an arrow with a solid, dark arrowhead denotes a path 
segment, and the arrow may be heavy to mark a level segment. In the figures 
where the direction of the edges are important, we again use arrows to mark 
the upward direction, but we make the arrowheads V-shaped ( "open" ) in this 
case to differentiate the edges from the segments. Dotted lines are used to 
show level lines in a face. 

A path P from s to t on the terrain is descending if the ^-coordinate of a point 
p never increases while we move p along the path from s to t. We assume that 
all paths and segments in our discussion are directed. Our discussion relies on 
the following known [1,17] properties of an SDP: 
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-> Edge with upward direction 
Level line 

Free segment 
Constrained segment 



Fig. 1. General legend for the figures in this paper 
Lemma 1 Any subpath of an SDP is an SDP. 

Lemma 2 An unfolded SDP is not always a straight line segment. 

Lemma 3 The intersection of an SDP P with a face of the terrain is either 
empty or a line segment. 

2.2 Related Work 

The SDP problem was introduced by de Berg and van Kreveld [0], who gave 
an algorithm to preprocess a terrain in O(nlogn) time so that it can be 
decided in O(logn) time if there exists a descending path between any pair 
of vertices. They did not consider the length of the path, and left open the 
problem of finding the shortest such path. Roy, Das and Nandy [17] solved 
the SDP problem for two special classes of terrains. For convex (or concave) 
terrains, they use the continuous Dijkstra approach to preprocess the terrain in 
O(n^logn) time and O(n^) space so that an SDP of size k can be determined 
in 0{k + logn) time. For a terrain consisting of edges parallel to one another, 
they find an SDP in O(nlogn) time by transforming selected faces of the 
terrain in a way that makes the unfolded SDP a straight line segment. Roy [f 6] 
has recently improved this running time to 0{n), by replacing a sorting step 
in the previous algorithm with a divide-and-conquer technique. Ahmed and 
Lubiw [1] examined the basic properties of SDPs that show the similarities 
and the dissimilarities between SDPs and shortest paths, and indicated why 
a shortest path algorithm like the continuous Dijkstra approach cannot be 
used directly to solve the SDP problem on general terrains. They also gave an 
0{n^'^\og{Y)) time algorithm that finds a (1 + e)-approximate SDP through 
a given sequence of faces. Their algorithm first formulates the problem as a 
convex optimization problem, which is then solved using a standard technique 
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in convex programming. In a more recent work Ahmed and Lubiw [3] gave a 
full characterization of the bend angles of an SDP, which shows that the bend 
angles along an SDP follows a generalized form of Snell's law of refraction 
of light. This result implies that computing an exact SDP is not easy even 
when we know the sequence of faces used by the SDP, due to numerical issues 
similar to the ones faced by Mitchell and Papadimitriou while computing a 
shortest path in the Weighted Region Problem [13, Section 8]. 

It was Papadimitriou [11] who first introduced the idea of discretizing space 
by adding Steiner points and approximating a shortest path through the space 
by a shortest path in the graph of Steiner points. He did this to find a shortest 
obstacle-avoiding path in 3D — a problem for which computing an exact solu- 
tion is NP-hard [7]. On polyhedral surfaces, the Steiner point approach has 
been used in approximation algorithms for many variants of the shortest path 
problem, particularly those in which the shortest path does not unfold to a 
straight line segment. One such variant is the Weighted Region Problem [13]. 
In this problem, a set of constant weights is used to model the difference in 
costs of travel in different regions on the surface, and the goal is to minimize 
the weighted length of a path. Mitchell and Papadimitriou [13] used the con- 
tinuous Dijkstra approach to get an approximate solution in O (n^log(-)) 



time. Following their result, several faster approximation schemes [4,5,6,8,22] 
have been devised, all using the Steiner point approach. The Steiner points 
are placed along the edges of the terrain, except that Aleksandrov et al. [6] 
place them along the bisectors of the face angles. A comparison between these 
algorithms can be found in Aleksandrov et al. [U]. 

One generalization of the Weighted Region Problem is finding a shortest 
anisotropic path [15], where the weight assigned to a region depends on the 
direction of travel. The weights in this problem capture, for example, the effect 
the gravity and friction on a vehicle moving on a slope. Lanthier et al. [12], 
Sun and Reif [21] and Sun and Bu [19] solved this problem by placing Steiner 
points along the edges. 

Note that all the above-mentioned Steiner point approaches place the Steiner 
points in a face without considering the Steiner points in the neighboring 
faces. This strategy works because we can travel between any two points in 
a face. In the case of shortest anisotropic paths, the straight-line path may 
be in a forbidden direction, but it is almost always assumed that the allowed 
directions permit a zigzag path to any destination (like tacking against the 
wind in a sailboat). The one exception is that Sun and Reif [21] consider the 
Anisotropic Path Problem where a set of non- adjacent faces have directions 
that are unreachable even with zigzagging. Their solution involves propagat- 
ing extra Steiner points across each of these partially traversable faces. For 
the SDP problem, ascending directions are unreachable in every face, which 
necessitates our non-local strategy of placing Steiner points. 
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2. 3 The Bushwhack Algorithm 



To compute a shortest path in the graph of Steiner points in a terrain we 
use a variant of Dijkstra's algorithm developed by Sun and Reif [20]. Their 
algorithm, called the Bushwhack algorithm, achieves 0(1^1 log |y|) running 
time by utilizing certain geometric properties of the paths in such a graph. 
The algorithm has been used in shortest path algorithms for the Weighted 
Region Problem [6,22] and the Shortest Anisotropic Path Problem [21]. 





(a) 



(b) 




Fig. 2. Maintaining the list Ie,e' in the Bushwhack algorithm 

The Bushwhack algorithm relies on a simple, yet important, property of short- 
est paths on terrains: two shortest paths through different face sequences do 
not intersect each other at an interior point of a face. As a result, for any two 
consecutive Steiner points Ui and U2 on edge e for which the distances from 
s are already known, the corresponding sets of "possible next nodes on the 
path" are disjoint, as shown using shading in Figure 2(a). This property makes 
it possible to consider only a subset of links at a Steiner point v when expand- 
ing the shortest path tree onwards from v using Dijkstra's algorithm. More 
precisely. Sun and Reif maintain a dynamic list of intervals /e,e' for every pair 
of edges e and e' of a common face. Each point in an interval is reachable from 
s using a shortest path through a common sequence of intermediate points. 
For every Steiner point v in e with known distance from s, /e,e' contains an 
interval of Steiner points on e' that are likely to become the next node in the 
path from s through v. The intervals in /e.e' are ordered in accordance with 
the ordering of the Steiner points v on e, which enables easy insertion of the 
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interval for a Steiner point on e whose distance from s is yet unknown. For 
example, right after the distance of M4 from s becomes known (i.e., right after 
U4 gets dequeued in Dijkstra's algorithm) as shown in Figure 2(b), the inter- 
val of the Steiner points on e' that are closer to than to any other Steiner 
points on e with known distances from s can be computed in time logarithmic 
in the number of Steiner points on e', using binary searches (Figure 2(c)). 
Let X denote this interval for ease of discussion. The Bushwhack algorithm 
considers only the Steiner points lying in interval X as the possible next nodes 
on the path to U4, while Dijkstra's algorithm tries all the Steiner points on 
e'. Another difference between these two algorithms is that after gets de- 
queued, Dijkstra's algorithm enqueues each Steiner point (or sifts it upward 
in the queue if it was already there) of X. But in the Bushwhack algorithm, 
only the Steiner point u'^ E X that is nearest from U4 is enqueued (or sifted 
upward); other Steiner points in X are considered later on if necessary, one by 
one and in order of their distances from u'^. Since u'^ can be located in interval 
X in constant time, each iteration of the Bushwhack algorithm takes 0(|y|) 
time, resulting in a total running time of 0(|1^| log |^|). 



2.4 Placing the Steiner Points 



Our approximation algorithms work by first discretizing the terrain with many 
Steiner points along the edges, and then determining a shortest path in a di- 
rected graph in which each link connects a pair of vertices or Steiner points in 
a face of the terrain in the descending (more accurately, in the non-ascending) 
direction. Although the idea is similar to other Steiner point approaches dis- 
cussed in Section 2.2, there are two aspects of the SDP problem that make 
our approach quite different from previous Steiner point approaches. 




Fig. 3. Problems with independently-placed Steiner points. 

First, because of the nature of the SDP problem, we have to position the 
Steiner points quite differently from the Steiner point approaches discussed in 
Section 2.2. In particular, we cannot place Steiner points in an edge without 
considering the heights of the Steiner points in other edges. More elaborately, 
for each Steiner point p in an edge, if there is no Steiner point with height 
h{p) in other edges of the neighboring faces, it is possible that a descending 
path from s to v through Steiner points does not exist, or is arbitrarily longer 
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than the SDP. For example, consider the SDP P = {s,pi,p2,P3, v) in Figure 3, 
where for each i G [1, 3], Qi, q[ and q'l are three consecutive Steiner points with 
h{qi) > h{q[) > h{q") such that is the nearest Steiner point above Pi. Note 
that in this figure the faces have been unfolded onto a plane, and that pi 
and q'l are the same point. There is no descending path from s to ^; through 
the Steiner points: we must cross the first edge at q'l or lower, then cross the 
second edge at or lower, and cross the third edge at q'^ or lower, which puts 
us at a height below h{v). Another important observation is that even if a 
descending path exists, it may not be a good approximation of P. In Figure 3, 
for example, if we want to reach instead a point v' slightly below v, P' would 
be a feasible path, but the last intermediate nodes of P and P' are not very 
close. We can easily extend this example to an SDP P going through many 
edges such that the "nearest" descending path P' gets further away from P at 
each step, and at one point, P' starts following a completely different sequence 
of edges. Clearly, we cannot ensure a good approximation by just making the 
Steiner points on an edge close to each other. 

To guarantee the existence of a descending path through Steiner points that 
approximates an SDP from s to any vertex, we have to be able to go through 
the Steiner points in a sequence of faces without "losing height", i.e., along 
a level path. We achieve this by slicing the terrain with a set of horizontal 
planes, and then putting Steiner points where the planes intersect the edges. 
The set of horizontal planes includes one plane through each vertex of the 
terrain, and other planes in between them that are close enough to guarantee 
a good approximation ratio. Our two algorithms, discussed in Sections 3 and 4, 
differ from each other in the manner the positions of the horizontal planes are 
determined. 

The second issue is that the previous Steiner point approaches relied on the 
property that shortest paths in the Weighted Region Problem or in the Short- 
est Anisotropic Path Problem cannot become very close to a particular vertex 
more than once. This property does not hold for shortest paths in the SDP 
problem. In fact, it is possible to construct a terrain where an SDP becomes 
very close to a vertex v as many as 0{n) times, moving far away from v after 
every visit of the vicinity of v. Consider the terrain in Figure 4(a) which con- 
sists of the triangular faces of a pyramid with a star-shaped base. The points 
s and t have the same height, so the SDP P from s to t must consist of level 
segments. Moreover, P consists of 0{n) segments in the figure. Figure 4(b) 
shows the faces used by P after unfolding them onto a plane. By moving the 
convex vertices at the base away from the "center" of the base while keep- 
ing them on the same plane, we can make the points of P that are far away 
from V move even further away from v. Clearly it is possible to make P enter 
and leave a region close to v as many as 0{n) number of times. Because of 
such a possibility with an SDP, the analysis of our Steiner point approach is 
completely different from previous approaches. 
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(b) 

Fig. 4. An SDP that comes close to a vertex 0{n) number of times 
3 Discretizing using Uniform Steiner Points 

In our first algorithm the Steiner points on each edge are evenly spaced. To 
determine their positions, we first take a set of horizontal planes such that 
any two consecutive planes are within distance 6 of each other, where 6 is 
a small constant that depends on the approximation factor. We then put a 
Steiner point at the intersection point of each of these planes with each of 
the terrain edges. One important observation is that this scheme makes the 
distance between consecutive Steiner points on an edge dependent on the 
slope of that edge. For instance, the distance between consecutive Steiner 
points is more for an almost-level edge than for an almost vertical edge. Since 
6 is the largest acute angle between a non-level edge and a vertical line, it 
can be shown that the distance between consecutive Steiner points on a non- 
level edge is at most 6 sec 9 (Lemma 5). Because of the situation depicted 
in Figure 3, we cannot place extra Steiner points only on the edges that 
are almost level. We guarantee a good approximation ratio by choosing 6 
appropriately. More precisely, we make sure that d sec 9 is small enough for 
the desired approximation ratio. Note that we can put Steiner points on a 
level edge without considering heights, since a level edge can never result in 
the situation depicted in Figure 3 (because all the points in such an edge have 
the same height). 
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3.1 Algorithm 



Our algorithm runs in two phases. In the preprocessing phase, we place the 
Steiner points, and then construct a shortest path tree in the corresponding 
graph. During the query phase, the shortest path tree gives an approximate 
SDP in a straightforward manner. 



3.1.1 Preprocessing Phase 

Let 6 = '^'^^^ ^ . We subdivide every non- level edge e of the terrain by putting 
Steiner points at the points where e intersects each of the following planes: 
z = j6 for all positive integers j, and z = h{x) for all vertices x of the terrain. 
We subdivide every level edge e by putting enough Steiner points so that the 
length of each part of e is at most 5 sec 6. Let V be the set of all the vertices 
and Steiner points in the terrain. We then construct a weighted directed graph 
G = (y, E) as follows, starting with i? = 0. For every pair (x, y) of points in 
V adjacent to a face / of the terrain, we add to £^ a directed link from x 
to y if and only if h{x) > h{y) and xy is either an edge of the terrain or a 
segment through the interior of /. Note that we do not add a link between 
two points on the same edge unless both of them are vertices. Each link in E 
is assigned a weight equal to the length of the corresponding line segment in 
the terrain. Finally we construct a shortest path tree T rooted at s in G using 
the Bushwhack algorithm. 

Note that we are mentioning set E only to make the discussion easy. In prac- 
tice, we do not construct E explicitly because the neighbors of a node x E V in 
the graph are determined during the execution of the Bushwhack algorithm. 




10 



When the query point v is a. node of G, we return the path from s to in T as 
an approximate SDP. Otherwise, we find the node u among those in V lying 
in the face(s) containing v such that h{u) > h{v), and the sum of the length 
of the path from s to w in T and the length of the segment uv is minimum. 
We return the corresponding path from s to as an approximate SDP in this 
case. To elaborate more on the latter case, let U be the set consisting of the 
nodes u E V with the following properties: 

(i) u and v lie in a common face, and 

(ii) h{u) > h{v). 

It is easy to see that if v is an interior point of a face, then all the nodes in 
U lie on at most three edges of that face (Figure 5(a)). Otherwise, v is an 
interior point of an edge, and there are at most four edges on which the nodes 
in U can lie (Figure 5(b)). Since we already know the length of an SDP from 
s to any u E U, we can find in \U\ iterations the node u E U that minimizes 
the length of the path constructed by concatenating the segment uv at the 
end of the path from s to m in T. The corresponding path is returned as an 
approximate SDP. 



3.2 Correctness and Analysis 

For the proof of correctness, it is sufficient to show that an SDP P from s to 
any point v in the terrain is approximated by a descending path P' such that 
all the segments of P', except possibly the last one, exist in G. We show this 
by constructing a path P' from P in the following way. Note that P' might 
not be the path returned by our algorithm, but it provides an upper bound 
on the length of the returned path. 

Let P = {s = po,pi,p2, . . . ,pk,v = Pk+i) be an SDP from s to v such that pi 
and pi^i are two different boundary points of a common face for all i G [0, A;— 1] , 
and pk and pk+i are two points of a common face. For ease of discussion, let Ci 
be an edge of the terrain through pi for all i G [l,k] [ci can be any edge through 
Pi if Pi is a vertex). Intuitively, we construct P' by moving each intermediate 
node of P upward to the nearest Steiner point. More precisely, we define a 
path P' = {s = Pq,p'i,P2, ■ ■ ■ ,p'k^v = p'k^i) as follows. For each i e [1, k], let 
= Pi if Pi is a vertex of the terrain. Otherwise, let p'^ be the nearest point 
from Pi in V n Ci such that h{p'i) > h{pi). Such a point always exists in V 
because pi is an interior point of Cj in this case, and it has two neighbors x 
and y in V n Ci such that h{x) > h{pi) > h{y). Note that each node of P' 
except possibly the last one is either a vertex or a Steiner point. 

Lemma 4 Path P' is descending, and the part of P' from s to p'^, exists in G. 
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PROOF. We prove that P' is descending by showing that h{p'j) > ^(Pi+i) 
for every i G [0,k]. We have: h{p'^) > h{pi+i), because h{p'j) > h{pi) by the 
definition of p[, and h{pi) > as P is descending. Now consider the 

following two cases: 

Case 1: p'^^^ = pi+i or Cj+i is a level edge. In this case, /^(Pi+i) = h{pi+i). It 
follows from the inequality h{p'j) > h{pi+i) that h{p[) > h{p[_^_^). 

Case 2: p'-^^ ^ pi^i and ej+i is a non-level edge. In this case, there is either 
one or no point in Cj+i at any particular height. Let p'-_^i be the point in 
ej+i such that /i(jo^Yi) — ^ip'i)^ if ^'^ such point exists, let p'-_^-^ be the 
upper vertex of Cj+i. In the latter case, we can infer from the inequality 
h{p'i) > h{pi^i) that h{p[) > h{p'l_^^. Therefore we have h{p'^) > h{p'-^^) 
in both cases. Since p'-^^ E V H Cj+i, the definition of p[_^_-^ implies that 
Hi/U) > HpU,). So, hip'^ > h{p[^,). 

Therefore, P' is a descending path. 

To show that the part of P' from s to p'f. exists in G, it is sufficient to prove 
that p'iP'ij^i G E for all i G [0, A; — 1], because both p[ and p'^j^^ are in V by 
definition. We have already proved that h{p'j) > Since p[ and are 

boundary points of a common face by definition, p'iP'i^i ^ E only in the case 
that both of and lie on a common edge, and at most one of them is 
a vertex. We show as follows that this is impossible. When both pi and pi+i 
are vertices of the terrain, both p'^ and p'^_^_i are vertices. When at least one 
of Pi and Pi+i is an interior point of an edge, they cannot lie on a common 
edge [1, Lemma 3]; therefore, both of and Pi_^_i cannot lie on a common edge 
unless both of p[ and p[_^_^ are vertices. So, this is impossible that both p[ and 
p^_^i lie on a common edge, and at most one of them is a vertex. Therefore, 

p'iP'i+^eE. □ 

Lemma 5 For all i E [l,k], \piPi\ < SsecO. 



PROOF. 




Fig. 6. Bounding \pip[\ when pi ^ p'- and Cj is a non-level edge 
When Pi = p'j^, \piPi\ = < 5 sec 6*. When pi ^ p[, and is a level edge. 
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\piP'i I < 5 sec 6 by construction. We will now focus on the case pi ^ p\ and Cj 
is a non-level edge. 

Consider the vertical plane containing the edge Cj. Construct a line vertically 
upward from pj to the point where = h{j^^ (Figure 6). Let Qi be the 
angle Aqi'p(p\. Since h{(ii) = h{p'^ > h{pi), 9i is an acute angle, and hence 
9 > 9i, which implies: 

cos ^< cos = 1^*^ =^ \piPi\ < iQiPil sec0 . 
\PiPi\ 

As QiPi is a vertical line, 

\qiPi\ = h{qi) - h{pi) = h{pi)-h{pi) < 5 
by construction, and therefore, \pip'i\ < \qiPi\ sec 9 < 6 sec 9. □ 
Lemma 6 Path P' is a {1 + e)- approximation of P. 



PROOF. When k = implying that P does not cross an edge of the terrain, 
we have P = (s, v) = P' which proves the lemma trivially. We will now focus 
on the case k > Q. 

The length of P' is equal to: 



E \p'iP'i+i\ < E [\p'iP^\ + \PiPi+i\ + \p 

i+ip'ij^i I ] (from triangle inequality) 

j=0 i=0 

k k 



= E \PiPi+i I + 2 E \PiPi I (since po = p'q and pk+i = p'k+i) 

1=0 i=l 
k k 

< ^ \PiPi+i I + 2 E sec 6^ (Lemma 5) 

1=0 i=l 
k 

< E bi^'^+i I + 2^'*^ sec 6* . 



i=0 



Because the number of faces in the terrain is at most 2n, and P has at most 
one segment in each face (Lemma 3), we have: k < 2n. Therefore, 



k k 

E \p'iP'i+i I < E \PiPi+i I + 4?^^ sec 9 

i=0 i=0 
k 

= J2\PiPi+'i- \ + ' 
i=0 
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from the definition of S. Because k > 0, pi lies on the edge opposite to po in 
the face containing both po and pi, and therefore, h < \poPi\ < Y.i=o \PiPi+i\- 
So, 

k k 

H \p'iP'i+i\ < (1 + e) H \PiPi+i\ ■ 

i=Q i=0 



Since P' is descending (Lemma 4), it follows that P' is a (1 + e) -approximation 
of P. □ 



Lemma 7 Let X = (^^^ sec 9. Graph G has less than i^^^ nodes and O 
links. Moreover, it has less than nodes along any edge of the terrain. 



PROOF. We will first prove the last part of the lemma. For each edge e of 
the terrain, the number of Steiner points corresponding to the planes z = j6 
is at most f — 1, and the number of Steiner points corresponding to the planes 
z = h{x) is at most n — 2. So, 

1^ n e| < (^j - 1^ + {n-2) + 2 < j + n = in (^^^ (^-^ sec 6* + n, 



because 6 = Since (f) (A sec^ > 1, we have: 



y n e < 5n — - sec6' = 



hj \e 



We will now compute \V\ and \E\. Let c = for ease of discussion. Using 
the fact that the number of edges is at most 3n, we have: 

V < 3nc = . 



For each face / of the terrain, there are less than 3c points in V f] f, and 
each such point has less than 2c neighbors in / (more precisely, in the induced 
subgraph G[V fl /]). So, the number of directed links in E contributed by / 
is less than 6c^, and this bound is tight for a level face. Because there are at 
most 2n faces, 

1^1 < Unc" = O ( . □ 
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Theorem 8 Let X = ^-^^ sec 6. Given a vertex s, and a constant e G (0, 1], we 
can discretize the terrain with -^^'"^"^ Steiner points so that after a preprocessing 
phase that takes O (j^^ log (~)) time for a given vertex s, we can determine 
a (1 + e) - approximate SDP from s to any point v in: 

(i) 0{n) time if v is a vertex of the terrain or a Steiner point, and 

(ii) O (^^^ time otherwise. 



PROOF. We first show that the path P" returned by our algorithm is a 
(1 + e)-approximation of P. Path P" is descending because any path in G is 
a descending path in the terrain, and the last segment of P" is descending. It 
follows from the construction of P" that the length of P" is at most that of 
P', and hence by Lemma 6, P" is a (1 + e)-approximation of P. 

As we have mentioned before, we do not construct E explicitly because the 
neighbors of a node x & V in the graph are determined during the execution of 
the Bushwhack algorithm. As a result, the (implicit) construction of G takes 
0(|y|) time. It follows from the running time of the Bushwhack algorithm 
(discussed in Section 2.3) that the preprocessing time of our algorithm is: 

0{\V\\og\V\) = of^log^^^^ 



by Lemma 7. 

During the query phase, if v is a vertex of the terrain or a Steiner point, the 
approximate path is in the tree T. Because the tree has height 0{n), it takes 
0{n) time to trace the path. Otherwise, v is an interior point of a face or an 
edge of the terrain. The last intermediate node u on the path to v is a vertex 
or a Steiner point that lies on the boundary of a face containing v. If v is 
interior to a face [an edge] , there are 3 [respectively 4] edges of the terrain on 
which u can lie. Thus there are O choices for u by Lemma 7, and we try 
all of them to find the best approximate path, which takes: 

O i^] + Gin) = G 



time. □ 



Note that the space requirement of our algorithm is 0(| V|) = G j since we 
are not storing E explicitly. Also note that using Dijkstra's algorithm with a 
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Fibonacci heap [11] instead of the Bushwhack algorithm yields an even simpler 
algorithm with a preprocessing time of 0(1^1 log \V\ + \E\) = O (n^ f^) ) • 



4 Discretizing using Steiner Points in Geometric Progression 

Unlike our first algorithm where the Steiner points on each edge are evenly 
spaced, our second algorithm places them non-uniformly along the edges. The 
Steiner points we use here are of two kinds. We first place Steiner points in 
"geometric progression" along the edges, as done by Aleksandrov et al. [1]. We 
call these points primary Steiner points. Then we place more Steiner points, 
called isohypse Steiner points, to guarantee that for every descending path in 
the terrain there exists a descending path through the Steiner points. Although 
the number of Steiner points used in this technique is more than in our first 
algorithm, the running time of the resulting algorithm no longer depends on 
the slope of the edges. 

4-1 Algorithm 

4.1.1 Preprocessing Phase 

The primary Steiner points are placed in such a way that for each vertex v of 
an edge e, there is a set of primary Steiner points whose distances from v form 
a geometric progression. Although the distance between a pair of consecutive 
Steiner points on e increases as we move away from f , we can still guarantee a 
good approximation ratio. This is because intuitively the length of a segment 
connecting two edges adjacent to v increases as we move the segment away 
from V — see Lemma 11 for a more precise statement. One observation is that if 
we want to maintain the geometric progression of the distances for the Steiner 
points very close to v, we would need infinitely many Steiner points near v. 
To avoid this problem, we do not put any primary Steiner points in a small 
region near v. 

Before going into further details, we will define a few constants for ease of 
discussion. Let ^1 = ^1 and <^2 = |x- The constant 5i will define a region near 
V where we do not put any primary Steiner points, while 82 will determine the 
distances between consecutive primary Steiner points outside that region. 

Definition 9 (Vicinity of a Vertex) In a face f incident to a vertex v, 
let pi and p2 be two points lying on two different edges of f at v such that 
\vpi\ = \vp2\ = Si. Clearly, AvpiP2 is an isosceles triangle. The vicinity of v 
is defined to be the union of all such isosceles triangles around v (Figure 7). 
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Fig. 7. Vicinity of a vertex 

Note that the vicinities of any two vertices Vi and V2 are mutually disjoint 
because 5i < | < ^i^^. 

In the preprocessing phase, we determine the positions of the Steiner points 
as follows. First, on every edge e = V1V2 we place primary Steiner points at 
points p E e such that \pq\ = Si{l + 62)^ for q G {vi, V2} and i G {0, 1, 2, . . .}. 
Then we add up to 3n isohypse Steiner points for each primary Steiner point 
and for each vertex, as follows. For every non-level edge e, and every point p 
that is either a primary Steiner point or a vertex, we place an isohypse Steiner 
point at the point where e intersects the horizontal plane through p (i.e., the 
plane z = h{j>)). 

After placing the Steiner points, we construct a weighted directed graph G = 
(V, E) and then construct a shortest path tree T rooted at s in G in the same 
way as in our first algorithm (Section 3.1.1). 

4- 1-2 Query Phase 

The queries are handled in exactly the same manner as in Section 3.1.2. 
4-2 Correctness and Analysis 

For the proof of correctness, we follow the same approach used in Section 3.2: 
given an SDP P, we first construct a path P' by moving each intermediate 
node of P upward to the nearest Steiner point, and then show that P' is 
descending and that it approximates P. This proves the correctness of our 
algorithm because the path returned by our algorithm is not longer than P'. 

Let P = {s = Po,Pi,P2, ■ ■ ■ ,Pk,v = Pk+i) be an SDP from s to v such that Pi 
and are two different boundary points of a common face for all i G [0, A;— 1] , 
and pk and pk+i are two points of a common face. Let Cj be an edge of the 
terrain through pi for all i G [1, A;]; Cj can be any edge through pi if pi is a 
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vertex. Now define path P' = {s = Pq,p'i,P2, ■ ■ ■ ,p'kiV = p'^^i) as follows: for 
each i G [1, A;], let = Pi if pt is a vertex of the terrain; otherwise, let be 
the nearest point from pi in V f] Cj such that h{p'j) > h{pi). 

Lemma 10 Path P' is descending, and the part of P' from s to p'^ exists in 
G. 

PROOF. The proof is exactly the same as in Lemma 4. □ 

Lemma 11 For all i G [l,k] such that pi is not inside a vertex vicinity, 

\PiPi\ < a \Pi-^Pi\ ■ 




(a) (b) 
Fig. 8. Bounding \pip'j\ when the face angle at vi is (a) acute and (b) obtuse 

PROOF. If Pi coincides with p^, the lemma follows trivially as \pip'i\ = 0. We 
will now focus on the case when these two points do not coincide. Since Pi is 
not inside a vertex vicinity, there is another Steiner point p'{ in Cj such that 
and p'- lie on the opposite sides of Pi. Let Vi be the common vertex of ej_i 
and Ci, Wi be the other vertex of Cj, and qi and w'i be two points in ej_i such 
that piQi -L Ci-i and Wiw'i ± ej_i. Figure 8 depicts these vertices and points, 
for both the cases that the face angle at Vi is (a) acute and (b) obtuse. 

We will first show that \PiPi\ < S2\viPi\, and then prove the lemma using a 
property of similar triangles. We have two cases as follows. If \viPi\ < \viPi\, 
then by construction: 



M < il + S2)\Vip'^\ 

^ hp'il - \viPi\ < S2 \viPi\ 

^ \p'iPi\ < S2 \ViPi\ < 82 \ViPi\ , 

since Pi lies strictly in between p[ and p'l- On the other hand, if \viPi \ > \vip'i 
then by construction: 
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\Vip'l\ < {l + 52)\v^p[\ 
-^\vip'l\ - \Vip[\ < 62 \Vip[\ 

=^\p'iPi\ < S2\v,p[\ < 62\ViPi\ , 

since pi lies strictly in between p[ and p'-. In both cases, \p[Pi\ < Si2\viPi 
We have: 

\PiP'i\ < \PiPi\ < ^2 \ViPi\ 
\ViPi\ 



h \qiPi 



iQiPi 



■^2 \QiPi\ ■ I \ \ (since AviPiQi and AviWiw'^ are similar) 



<^2 \qiPi\ ■ T 

h 

eh L 

■ iQiPil ■ ~r (from the definition of ^2) 



<^\pi-iPi\ {since \pi-ipi\ > \qiPi\) . □ 



6L '^ ^ ' h 
e 

6 

Lemma 12 For all z G [i-,k] such that pi is on or inside a vertex vicinity, 
, eh 

PROOF. If Pi is a vertex, the lemma follows trivially since p'^ = pi in this 
case. If Pi is not a vertex, let be the edge containing pi, and Vi be the vertex 
whose vicinity contains Pi. It is not hard to see that Vi is a vertex of Cj because 
^1 is strictly less than h. Let be the primary Steiner point on Cj which lies 
at distance 61 from Vi. Clearly pi lies in line segment ViQi. Now p'- cannot be 
outside line segment ViQi because otherwise we would have chosen either Vi or 
Qi as p'i- As a result, Pi also lies in hne segment ViQi. Therefore, 

I /I I I e 
iPiPil < \Viqi\ = Oi = — . □ 
on 

Lemma 13 Path P' is a [1 + e) - approximation of P. 
PROOF. The length of P' is equal to: 



Yl \p'iP'i+i\ < H {\p'iPi\ + \PiPi+i\ + \p i+iP'i+i I ) (from triangle inequality) 

i=0 1=0 
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= E \PiPi+l 
i=0 

k 

< E \PiPi+l 
i=0 

k 

= E \PiPi+l 
i=0 

k 

< E \PiPi+l 

k 

< J2 \PiPi+l 
1=0 



2 J2 \PiPt I (since Pq = p'o and pk+i = p'k+i) 



i=l 



2^ 7: \Pi~iPi\ + 7^ (by Lemmas 11 and 12) 



i;Y.\Pi-^Pi\ 
j=l 



e/i 
6n 

ehk 
3n 



3n 
2e/i 



since k < 2n because the number of faces in the terrain is at most 2n, and P 
has at most one segment in each face (Lemma 3). Assuming that P crosses at 
least one edge of the terrain (otherwise, P' = {s,v) = P), Yli=o \PiPi+i\ ^ 
and therefore: 



k k / e\ 2e ^ 

J2\PiPi+l\<Y.\P^Pi+^\('^+ n) + -irY^lPiPi+^l 
i=0 i=0 \ -3/ 6 

. e 2e\^, 
1 + 3 + yj E \PiPi+i\ 

k 

= (l + e) Y.\PiPi+i\ ■ 

Because P' is descending (Lemma 10), it follows that P' is a (H-e)-approximation 
of P. □ 



Observation 14 For any real number x G (0, 1], log(l + x) > 



PROOF. 



log(l + x) =loge ■ logg(l + x) 

( x^ x^ x^ x^ x^ \ 

= log e \ X 1 1 h... 

^ \ 2 3 4 5 6 / 

= log e(x(l--^+ (- - + x^ (- - - 
^ \ \ 2j \3 4j V5 6 

Since each term in the outer parentheses of the last expression is strictly 
positive, we have: 
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log(l + a;)>loge(^x(^l-|)) > log e (^x (^1 - = . □ 

Lemma 15 Graph G has less than log (^) nodes and O log^ (^)) 

links. Moreover, it has less than log nodes along any edge of the ter- 

rain. 



PROOF. We will first compute an upper bound on the number of primary 
Steiner points, which will then be used to prove the lemma. 

Let He be the number of primary Steiner points on edge e. It is straightforward 
to see that Ue is at most 2j, where j is the largest integer satisfying the 
following inequality: 



^1(1 + 52^ < L 



l0g(l+52) 



Therefore, 



ne<'2j < 



log(l + 52) 
< 1^ log ^ (Lemma 14) 



■^2 log e \J), 
2 

41og,2-6L, fL-6n 
log 



eh \ eh 

16.64L, /6nL 



< ; — log 



h \ eh 

Since there are at most 3n edges in the terrain, the total number of primary 
Steiner points is at most SnUe, which is less than: 

50nL fQnL 
log 



eh \ eh 



We will now prove the last part of the lemma. For each point p that is either 
a primary Steiner points and a vertex, there is at most one node inV P^e for 
any edge e. This is obvious when p lies on e. On the other hand, if p does not 
lie on e, there is at most one isohypse Steiner point on e that corresponds to 
p. Using the above bound on the number of primary Steiner points, we have: 
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, 50riL, /6nL\ 
y n e < — ; — log — — + n 
eh \ eh 



<^log(^j (smcef >1. 



We will now compute \ V\ and \E\. Let c = -^^^ log for ease of discussion. 

Using the fact that the number of edges is at most 3n, we have: 

iT/i Q ISSn^L, fQnL\ 

\V\ < 3nc = log — — 

eh \ eh J 



Using the same argument we used in the proof of Lemma 7, we can say that 
the number of directed links in E contributed by each / of the terrain is less 
than 6c^. Because there are at most 2n faces, 

Theorem 16 Given a vertex s, and a constant e G (0, 1], we can discretize 
the terrain with at most ^^^|-^log {^^^ Steiner points so that after a prepro- 
cessing phase that takes O (^^^log^ (^)) ^^"^^ ^ given vertex s, we can 
determine a (1 + e)- approximate S DP from s to any point v in: 

(i) 0{n) time if v is a vertex of the terrain or a Steiner point, and 
(a) 0(^^log(^)) time otherwise. 



PROOF. The proof is the same as in Theorem 16 except that we use Lem- 
mas 13 and 15 instead of Lemmas 6 and 7 respectively. 



As in the case of our first algorithm, we can use Dijkstra's algorithm with 
a Fibonacci heap [11] instead of the Bushwhack algorithm to have an even 
simpler algorithm with a preprocessing time of O {^t|t log^ (^) ) • 



5 Conclusion 



It may appear that the running time can be improved by using the technique 
by Aleksandrov et al. [6] who place Steiner points along the bisectors of the face 
angles. Although the technique improves all previous results on the Weighted 
Region Problem, it cannot be used for the SDP problem very easily. The main 
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problem is that it is not clear how to prove the existence of a feasible path 
that approximates an SDP. 

When query point v is neither a vertex of the terrain nor a Steiner point, the 
query phase can be made faster by using a point location data structure on 
each face. Note that the Voronoi diagram on each face consists of hyperbolic 
arcs. 
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